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Abstract 

We present metallicities and their variations with different parameters for 36 high- 
latitude fields covering Galactic longitudes 0° < I < 360°. The metallicities for 
relatively short vertical distances (z < 2.5 kpc) show systematic fluctuations with 
Galactic longitude, similar to those of the thick-disc scaleheight, which may be 
interpreted as indicating a common origin, viz., the flare effect of the disc (Bilir et 
al. 2007). This suggestion is supported by the metallicity variations which we find 
as functions of radial distance. 

The metallicity variation at larger vertical distances (6.5 < z < 9.5 kpc) is small 
but monotonic. Three different vertical metallicity gradients could be detected: 
d[M/H]/dz = -0.22(±0.03), d[M/H)jdz = -0.38(±0.06), and d[M/H]/dz = 
-0.08(±0.07) dex kpc" 1 for the intervals z < 3, 3 < z < 5, and 5 < z < 10 
kpc, respectively. Since our data cover the whole Galactic longitude interval, the 
resulting metallicity gradients can be interpreted as properties of the larger-scale 
Galaxy. The first gradient confirms the dissipational formation of the disc at short 
z-distances. The steeper gradient corresponds to the transition region between dif- 
ferent population components of the Galaxy, and finally, the lowest value provides 
an adequate description of the inner-halo metallicity gradient. 
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1 Introduction 



The metallicity distribution of G-type stars can be used for the interpretation 
of the formation and chemical evolution of the related Galactic components. 
For example, a metallicity gradient for a particular component of our Galaxy 
may indicate that this component formed by dissipational collapse, while the 
absence of a metallicity gradient may provide strong evidence of a merger 
and/or accretion origin from - possibly numerous - fragments, such as dwarf- 
type gala xies. The pioneering work on the dissip ational collapse of the Galaxy 
is that of lEggen. Lynden-Bell fc Sandaga (1l962f ) who argued that the Galaxy 
collapsed in a free-fall time ( ~2x 10 8 yr). Now, w e know that the Galaxy col- 



(19851: iNorrisI ( 19861 



l apsed over ma n y Gyr (e.g. lYoshii &: Saio I (119791) ; iNorris. Bessell fc Pickles 



bandage fc Fouts I fll987h: ICarnev. Latham fc Laird 



( 119901 ) ; INorris fc Ryan I (Il99ll ) ; iBeers &: Sommer-Larsen I (119951 ) ) . On the other 



ha nd, the merger or accr etion scenario of galaxy formation was first advocated 
by ISearle fc Zinn I (119781 ) ; f or a survey of the more recent l i teratu re on the sub- 



ject, we refer the reader to iFreeman fc Bland-Hawthorn I (120021 ) . 



Extant data comprise vertical metallicity gradients in the range —0.40 < 
d[M / H]/dz < —0.20 dex kpc" 1 for relatively small distances from the Galac- 



tic plane, i.e., z < 4 k pc (ITrefzger et al. II1995I ; iKaraali et al. 112003 ; iDu et al 



2004 ; lAk et al. 1 120071 ). which support the picture of a dissipative collapse of 



the thin disc. For intermediate ^-distances, where the thick disc is dominant, 
the vertical metallicity gradient is low, d[M/H]/dz = —0.07 dex kpc -1 , and 
the radial gradient is only ma rginal, —0.02 < d[M/H]/dz < dex kpc~ ] 



( Rong. Buser fc Karaali 1120011). We quote also the works of IChiba fc Beers 



(2000) and Girard et al. ( 1996 ) where vertical rotational velocity gradients 
were cited. In the first work, the authors give A < > /A\Z\ = —52 ± 6 km 
s _1 kpc^ 1 for the halo stars whereas in the second work A < > /A|Z| = 
—30 d= 3 km s _1 kpc -1 , for the stars dominated by the thick disc red giants. 
One must keep in mind that the samples assumed for a specific population, 
such as thin and thick discs or halo, m a y be contaminated by other objects. 
The work of iGilmore. Wyse. fc Norris I (120021 ) where higher mean rotational 
velocities than expected were predicted is a good example. Contrary to the 
expectations of the authors, the data of stars a few kpc from the Galactic 
plane which were assumed as candidates of the thick disc showed the mean 
rotation velocity < V >~ 180 km s -1 instead of the < V >~ 100 km s -1 
one. The origin of most of the stars in the sample was a disrupted satellite. 
Another important point in treating the metallicity or velocity gradient of 
a specific population is the effect of the metal-week component of the thick 
disc. If the separation of stars into different populations would be based mainly 
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on their metallicities, the metal-week tail of the thick disc may cause some 
complications. 



In our recent papers, Galactic model parameters for 36 high-latitude fields 
have been estimated using photometric data from the Two Micron Sky Survey 
(2MASS) and from the Sloan Digital Sky Survey (SDSS). For the thin and 
the thick discs, useful data have been available in both photometric surveys 
( Cabrera-Lavers et al. 20071; Bilir et al. 2007 ); however, for the halo only 



SDSS data could be used ( Bilir et al. 20071 ). We determined the variations of 



both the thick-disc scaleheight and the axial ratio of the halo as functions of 
Galactic longitude, and we noticed that it is most likely this variation which 
gives rise to the rather widely differing numerical values obtained for these 
parameters by different researchers. 

Our work was based on de-reddened apparent u, g, r, i, z magnitudes, esti- 
mated absolute magnitudes, and distances for 2 164 680 stars in a total area 
of 831 deg 2 . In the present paper, we shall derive metal abundances for a sub- 
sample of these stars and combine them with the distances already estimated 
in our previous paper, in order to investigate the large-scale metallicity distri- 
bution, including vertical and radial metallicity gradients, and their possible 
dependence(s) on Galactic longitude. 



Ak et al. I fl2007h showed that in the anti-centre direction, the metallicity gra- 



dients corresponding to the transition region from the thick disc to the halo 
(3 < z < 5 kpc) are significantly different at the two different Galactic lati- 
tudes b = +45° and b = —45°, respectively. In the present study, the variation 
of metal- abundances and metallicity gradients with Galactic longitude, a trace 
of the triaxiality of the halo will emerge as a particularly interesting result. 

A brief description of the SDSS data and their reductions is given in Section 
2. In section 3 we present the metallicity calibration for 36 fields and the 
longitude variation of the Galactic metal abundance distribution, including 
the vertical and radial metallicity gradients. These results are discussed in the 
final Section 4. 



2 SDSS 



The SDSS is a large, international collaboration project set up to survey 10 
000 square-degrees of sky in five optical passbands and to obtain spectra of 
one million galaxies, 100 000 quasars, and tens of thousands of Galactic stars. 
The data are being taken with a dedicated 2.5-m telescope located at Apache 
Point Observatory (APO), New Mexico. The telescope has two instruments: 
a CCD camera with 30 2048x2048 CCDs in the focal plane and two 320 fiber 
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double spectrographs. The imaging data are tied to a network of brighter 
astrometric standards (which would be saturated in the main imaging data) 
through a set of 22 smaller CCDs in the focal plane of the imaging camera. 
An 0.5-m telescope at APO has been used to tie the imaging data to brighter 
photometric standards. 



The SDSS obtains images almost simultaneously in five broad bands (u, g, 
r, i and zj}\ centered at 3551, 4686, 6166, 7480 and 8932 A, respectively 
(IFukugita et al. Ill996l ). The imaging data are automatically processed through 
a series of software pipelines which find and measure objects and provide pho- 
tometric and astrometric calibrations to produce a catalogue of objects with 
calibrate d magnitudes, positi ons and structure information. The photometric 
pipeline ( Lupton et al. IlioOll ) detects the objects, matches the data from the 
five filters, and measures instrumental fluxes, positions, and shape parameters 
(which allows the classification of objects as "point source", -compatible with 
the point spread function-, or "extended"). The photometric calibration is ac- 
curate to roughly 2 per cent rms in the g, r and i bands, and 3 per cent in u and 
z, as determined by the co nstancy of stellar population colours (llvezic et al 



20041 ; iBlanton et al. II2005I ). while the a strometric calibr ation precision is bet- 



ter than 0.1 arcsec rms per coordinate (IPier et al.l 120031). The Data Release 5 
(DR5) imaging catalogue covers 8000 deg 2 (lAdelman-McCarthy et al. 1 120071 ) 
with a detection repeatability complete at a 95 per cent level for point sources 
brighter than the limiting apparent magnitudes of 22.0, 22.2, 22.2, 21.3 and 
20.5 for u, g, r, i and z, respectively. The data are saturated at about 14 mag 
in g, r and % and about 12 mag in u and z. 



2. 1 Data and reductions 



The data used in this work were taken from the SDSS (DR5) WEB serve 
for 36 high-latitude fields (60° < b < 65°) covering different Galactic longitude 
intervals throughout the full circle (0° < I < 360°). From the 2 1 64 680 stars 



with SDSS magnitudes u, g, r, i, and z used in our previous paper (IBilir et al 



20071). we identified G-type stars by transformation of the Neil Reid's WEB 



L. we ic 



page|_| {U—B)q and (B — V) colours for the G spectral type dwarfs to the (u — 
g)o and (g — r) colours for this work. The absolute magnitude range of these 
stars, 258 268 in number in a total area of 831 deg 2 , is 5 < M(g) < 6. The range 
of apparent g magnitudes is 15 < go < 22, as in our paper just cited. While 
the de-reddening of the apparent magnitudes, the determination of absolute 
magnitudes and the estimation of the distances relative to the Sun as well as 



1 Magnitudes in this paper are quoted in the ugriz system to differentiate them 
from the former ugriz system. 

2 http://www.sdss.org/dr5/access/index.html 

3 http:/ /www-int. stsci.edu/~inr/intrins. html 
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Fig. 1. Mean error versus go apparent magnitude for stars with six apparent mag- 
nitude intervals, (15, 16], (16, 17], (17, 18], (18, 19], (19, 20] and (20, 21]. The mean 
error for 15 < go < 18 where thick disc stars are dominant is less than ±0.1 dex. 



the ^-distances from the Galactic plane are explained in iBilir et al. I (120071 ). 



derivation of the metal-abundances is described in the following Section 3. 



3 Metallicity 



The m etallicities of the sample s tars w ere evaluated using the following equa- 
tion of iKaraali. Bilir &: Tuncel I (120051 ): 



[M/H] = 0.10 - 3.545 .43 - 39.63<^ 43 + 63.51<^ 43 . 



This equation was calibrated for the main sequence stars with 0.12 < (g—r) < 
0.95 which covers the (g — r)o colour indices of our sample. Here, #0.43 is the 
normalized UV-excess in SPSS photometry corresponding to £0.6 m the UBV 
photometry. IKaraali. Bilir &: Tuncel I (120051 ) give the range of the metallicity 
as —2.76 < [M/H] < 0.2 dex, corresponding to < <5 .43 < 0.33. Limitations 
of the photometric technique are unavoidable, especially for the metallicities 
of faint star. However, the thick disc stars occupy the intermediate apparent 
magnitude interval, i.e. 15 < go < 18, (IChen et al.ll200ll ) where the mean error 
for [M/H] is less than ±0.1 dex (Fig. 1). We emphasize that the parameters 
based on the determination of the metallicity are limited by the accuracy of 
metallicities estimates based on photometry alone, and that they should be 
supplemented with investigations of the soon to emerge spectroscopic deter- 
minations of parameters for the SDSS samples. 

The parameter range < #0.43 < 0.33 covers the metallicity interval —3 < 
[M/H] < 0.2 dex. 

The metallicity distribution for each field is given by a mean metal- abundance 
value determined for each of the following distance intervals (in kpc): (1.5,2], 
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Table 1 

The metallicity distribution as a function of distance from the Galactic plane, for 36 
high latitude (60° < b < 65°) fields at different Galactic longitudes (0° < I < 360°). 
Uncertainties (given in parentheses) refer to the last quoted digits. Distances in kpc. 



l/z 1.07 1.58 2.02 2.46 3.17 4.04 4.92 5.80 6.68 7.56 9.50 
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(2,2.5], (2.5,3], (3,4], (4,5], (5,6], (6,7], (7,8], (8,9], (9,12], (12,15]. These in- 
tervals and the corresponding centroid distances (r*) are the same as in 
Ak et al. ( 20071 ); however, the projection of r* onto the vertical direction, 
i.e. z = r*sin(6), is different due to the higher Galactic latitudes studied in 
this work. 



Results are presented in Table 1. Notice that the metallicity distributions in 
the distance intervals 3<r<4, 4<r<5, 5<r<6, and 6 < r < 7 kpc 
are rather flat. Hence, median values could be used in doing the statistics for 
these intervals, whereas modal means were preferred for the other distance 
intervals which showed a peak. 
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Fig. 2. Mean metallicity distribution as a function of Galactic longitude, for vertical 
distances 1 < z < 2.5 kpc. 

3.1 Metallicity variation with Galactic longitude 

The metallicity distributions for both (relatively) short and large z distances 
show systematic fluctuations. For example, the distribution for z = 2 kpc and 
the combined distribution for 1 < z < 2.5 kpc have maxima at intermediate 
longitudes, while minima exist at small and large longitudes (Fig. 2). These 
z distances correspond to the region where the discs - both the thin and 
thick discs - dominate the stellar distributions of the Galaxy. If compared to 
the typical error bars, the amplitude of the fluctuation is significant but only 
moderate, i.e. A[M/H] ~ 0.2 dex. However, if it is real, this feature should 
be related with the structure of the Galactic disc, probably with its warp and 
flare (see Section 4). For completeness, the longitude variation of Fig. 2 has 
been fitted by the following equation (solid curve): 

[M/H] = -0.67 - 0.06 sin[27r(0.0046Z + 0.3940)]. (2) 



The metallicity distribution for the large z distances is even more interesting. 
Combining the distributions for 6.5 < z < 9.5 kpc reveals a clear longitudinal 
metallicity distribution (Fig. 3). The metal abundance decreases monotoni- 
cally from -2.28 to -2.38 dex when one goes from the Galactic longitude 10° to 
340°. However, this metallicity distribution can also be considered flat within 
the errors. The interesting is that the cited longitudes, 10° and 340°, are close 
to the direction of the Galactic centre, where one expects the same metallicity. 
The distance interval 6.5 < z < 9.5 kpc covers the halo component. Hence, 
this unexpected finding should be explained by the structure of the halo. The 
equation fitted to the metallicity distribution as follows: 

[M/H] = -0.00024/ - 2.292. (3) 
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Fig. 3. Mean metallicity distribution as a function of Galactic longitude, for vertical 
distances 6.5 < z < 9.5 kpc. 

3.2 Vertical and radial metallicity gradients 

Figure 4 displays the histograms of the metallicities given in Table 1 for the 
36 fields in each of 11 z distance intervals covering the range 1 < z < 9.5 
kpc. Each histogram has been fitted by a Gaussian curve with modal mean 
adopted as equal to the mean metallicity of 36 fields at the corresponding z 
distance from the Galactic plane (Table 2). As expected, the mean metallicities 
decrease with increasing z distances, indicating a probable vertical metallicity 
gradient. We plotted the same data in Fig. 5 for investigating this process in 
detail. One can see three different trends in the figure: 1) for relatively short 
z distances, i.e. z < 3 kpc, the variation of [M/H] is rather smooth, 2) for 
intermediate z distances, i.e. 3 < z < 5 kpc, the variation is steeper but still 
smooth, and finally 3) for 5 < z < 10 kpc, the variation is flat, but fluctuations 
and error bars are larger. The metallicity gradient describing the first trend, 
d[M/H]/dz = —0.22 dex kpc -1 , is in agreement with the canonical metallicity 
gradients for the same z distances, and is the likely signature of this Galactic 
region's formation by a dissipative collapse. The description of the second 
trend in terms of a metallicity gradient, d[M/H]/dz = —0.38 dex kpc -1 , 
corresponds to the average metallicity difference between the two population 
components involved. Finally, the third trend, d[M/H]/dz = —0.08 dex kpc -1 , 
is the low metallicity gradient of the inner spheroid. 

We investigated the radial metallicity gradient for 36 fields as follows. It will be 
seen immediately that the radial metallicity gradient described and used in this 
work is different from the usual definition: although we defined metallicities 
for 36 fields for a specific z distance from the Galactic plane, the distances 
of the centers of these fields to the Galactic centre are, of course, different 
due to their different Galactic longitudes, 0° < / < 360°. The range of radial 
distances thus defined is 7 < R < 16 kpc. Metallicities for 36 fields in Table 
1 as functions of radial distance R are given in 11 panels in Fig. 6. From the 
definition, short radial distances correspond to fields in the general direction 
toward the centre of the Galaxy, whereas fields in the general direction toward 
the anti-centre have larger radial distances. 
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Fig. 4. Metallicity histograms for the mean metallicities of 36 fields. 
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Fig. 5. Mean metallicity for 36 fields as a function of vertical distance z. 

Different trends can be observed in different panels. The most conspicuous fea- 
tures in panel (a) are the following: the fields in both the centre and anti-centre 
directions are relatively richer in metallicity than the fields in other directions, 
and their corresponding error bars are larger; furthermore, the metallicity 
differences between the centre and anti-centre fields follow an axisymmetric 
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Table 2 

Mean metallicities for different distance intervals calculated from the metallicity 
distributions of 36 fields. Symbols: r: distance from the Sun; r*: centroid distance 
corresponding to the interval r\ — r2\ z: projection of r* onto the vertical direction; 
[M/H]: mean metallicity. 



r\ — T2 


* 

r 


z 


[M/H] 


(kpc) 


(kpc) 


(kpc) 


(dex) 


0.5-1.5 


1.20 


1.07 


-0.46(07) 


1.5-2.0 


1.78 


1.58 


-0.63(06) 


2.0-2.5 


2.28 


2.02 


-0.70(07) 


2.5-3.0 


2.77 


2.46 


-0.77(09) 


3.0-4.0 


3.57 


3.17 


-1.20(10) 


4.0-5.0 


4.56 


4.04 


-1.63(13) 


5.0-6.0 


5.54 


4.92 


-1.87(14) 


6.0-7.0 


6.54 


5.80 


-1.93(15) 


7.0-8.0 


7.53 


6.68 


-2.34(12) 


8.0-9.0 


8.53 


7.56 


-2.35(13) 


9.0-12.0 


10.71 


9.50 


-2.32(12) 



pattern. 

The relative metal overabundances in the anti-centre fields (i.e., for large radial 
distances) are even larger in panels (b) and (c). The metallicity variations in 
panels (d) and (e) are almost zero, and finally, the metallicity trend changes 
sign for higher z distances, i.e., 4 < z < 7 kpc (panels f-i). In short: fields in 
the centre direction of the Galaxy appear to be overabundant in metallicities 
relative to the other fields. Although this metallicity excess is relatively low, 
A[M/H] ~ 0.15 dex, it appears to be real because the present statistical 
analysis shows the metallicity trend behind it to be significant. 

For completeness, we should also mention that the last two panels, z = 7.56 
and z = 9.50 kpc, show zero variation in metallicity. 

We fitted linear equations (solid lines) to the metallicity distributions in Fig. 
6, and we plotted their slopes versus z distances in Fig. 7. Clearly, the metallic- 
ity variation exhibits different trends at different distances from the Galactic 
plane. The minimum at z ~ 5 kpc is the most conspicuous feature in this dia- 
gram, meaning that the slope of the radial metallicity gradient is significantly 
different at shorter and larger z distances. It thus seems that the Galactic 
components (thin and thick discs, and halo) do not have homogeneous struc- 
tures. 
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Fig. 6. Variation of the metallicities for 36 fields with radial distance R. 
4 Summary and Conlusion 

We evaluated metallicities for 36 high latitude (60° < b < 65°) fields with 
Galactic longitudes 0° < / < 360°, and we investigated their variations with 
different parameters. 

The metallicities are longitude dependent. For relatively short z distances from 
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Fig. 7. Radial metallicity gradient versus vertical distance z. 

the Galactic plane, z < 2.5 kpc, one can observe two minima at small and 
large Galactic longitudes (I = 80° and I = 300°), respectively, and a maximum 
at I = 190° (Fig. 2). For large distances, 6.5 < z < 9.5 kpc (Fig. 3), metallicity 
decreases monotonically with Galactic longitude. 



The metallicities of 36 fields were plotted as histograms, which were fitted 
by Gaussians giving a mean metallicity for each of 11 distinctive (representa- 
tive) values of z distance in the range 1 < z < 9.5 kpc. The resulting mean 
metallicity as a function of z reveals a trimodal behavior: (1) for z < 3 kpc, 
there is a metallicity gradient of d[M/H]/dz = —0.22 dex kpc -1 which is in 
agreement with the dissipational collapse picture of formation of the Galactic 
thin-disc component at relatively short z distances; (2) a steeper gradient, 
d[M / H]/dz = —0.38 dex kpc -1 , covering the distance interval 3 < z < 5 
kpc which, although dominated by the thick disc, probably accentuates the 
metallicity differences between the coexisting thin and thick discs and the 
halo, rather than it provides a pure fossil record of thick-disc formation; (3) 
the gradient d[M/H]/dz = —0.08 dex kpc -1 , for 5 < z < 10 kpc, is a typical 
value for the inner spheroidal component. 



These results, including the metallicity gradient d[M/ H]/dz = —0.26 dex 



kpc 1 


for z < 10 k 


pc, are all in fine agreement with the canonical litera- 


ture ( 


Trefzeer et al. 


1995: 


■lone:. Buser & Karaali 2001: Karaali et al. 


2003; 


Du et al. 


2004 




Ak et al. 


2007 


). By virtue of their systematic coverage of 



Galactic longitudes, they provide an important enhancement of our recon- 
naissance of Galactic structure which will be needed for an improved under- 
standing of the formation of the Milky Way Galaxy. 



Conversion of vertical (z) distances to radial (R) distances also allowed a study 
of the metallicity variations as functions of R over a range of almost 10 kpc. 
Again, clear (linear) systematic trends with different slopes can be identified 
behind the observed metall icity fluctuations, confirming results obtained by 



Allende Prieto et al. I (120061 1 from DR3 data for 22 700 F- and G-type stars. In 



particular, although these authors investigate a larger range in R, the metal- 
licity structures that they find for the thick disc and the halo (their figures 
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13) are almost the same as those presented here for z < 2.5 and 6.5 < z < 9.5 
kpc. 



Our re sults are reminiscent of the flare effect of the disc. According to lLopez-Corredoira et al 



( 120021 ). a flare occuring in the outer disc produces an increase in the scale- 
height as one moves radially outward. From their results we infer that the 
scaleheight of the thin disc sources is as large as 0.6 kpc for R = 10 kpc, so it 
is plausible that a mixture of sources from the thin and thick discs are found 
in the outer disc at the lowest height regimes plotted in Fig. 6 (panels a and 
b), producing higher metallicities than the expected from the contribution of 
the thick disc component alone. 

The longitudinal dependence of the metallicity can be then explained by this 
process, since different longitudes correspond to different radial distances for 
a specific distance from the Galactic centre. The anticentre direction corre- 
sponds to the larger distances from the Galactic centre, thus is in this longi- 
tude range where the higher metallicities must be obtained due to the effect 
of the flare, as it is shown in Fig. 2. Also the longitudes where either the 
maximum or the minimum of the metallicity distribution are observed coin- 
cide with the directions of max imum warp amplitude in the Galactic thin disc 



( lLopez-Corredoira et al. II2002I ). As this warp bends above and below the thin 
disc from its mean location, again mixing of stellar population from the thin 
and thick discs is expected at these longitudes: hence, observations of a varia- 
tion in the mean metallicity are not surprising. Finally, even minor metallicity 
differences between fields at longitudes I = 10° and I = 340° (both in the 
Galactic centre direction), they may result from a triaxial structure of the 
halo - or in other words: from the fact that the plane and the meridian of the 
disc are different from those of the halo. 
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